Cascade Boltzmann - Langevin approach to higher-order current correlations in 

diffusive metal contacts 
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The Boltzmann - Langevin approach is extended to calculations of third and fourth cumulants 
of current in diffusive-metal contacts. These cumulants result from indirect correlations between 
current fluctuations, which may be considered as "noise of noise". The calculated third cumulant 
coincides exactly with its quantum-mechanical value. The fourth cumulant tends to its quantum- 



mechanical value 
sign at eV ~ 20T. 



-e J/105 at high voltages and to a positive value 2e T/3R at V = changing its 
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The fluctuation - dissipation theorem relates the spec- 
tral density of noise in equilibrium with the correspond- 
ing linear response of the system. This is why the second 
cumulant of current presents new information about the 
system only under nonequilibrium conditions. Nonequi- 
librium noise in mesoscopic-systems has been extensively 
studied in the past decade.EI Unlike the second cumulant 
of current, its higher cumulants are not related with any 
average characteristics of the system even in equilibrium. 
Hence these quantities may be also of experimental in- 
terest. 

In recent years, higher-order correlations of currents 
in mesoscopic systems received considerable attention 
of theorists..-. This work was pioneered by Levitov 
and Lesovikp who calculated the distribution of charge 
transmitted through a single-channel quantum con- 
tact. Subsequently, these calculations were extended to 
phas&icoherent diffusive contacts in the zero-temperature 
limit u Very recently, the third cumulant of current was 
calculated for quasi-one-dimensional diffusive-jrnetal con- 
tacts for arbitrary temperatures and voltaeesu using the 
nonlinear cr-model and Keldysh formalismE 

Along with the fully quantum-mechanical approach, 
higher-order correlations have been studied using a 
semiclassical description, which ignores quantum in- 
terference of wave functions yet takes into account 
the Fermi statistics of electrons. Semiclassical calcu- 
lations of higher-order cumulants were performed for 
double-barrier structuresa using a master equation and 
for chaotic cavitiesQ based on the minimal correlation 
approach, a 

A semiclassical Boltzmann - Langevin approachta has 
been very successful in describing the second cumulant 
of current in diffusive contacts.£3 In particular, it gives 
the same value of the shot noise-|S7 = 2e//3 as a fully 
quantum- mechanical treatment O This is not surpris- 
ing because the corrections to semiclassical values from 
quantum-interference effects like weak localization and 
mesoscopic fluctuations ofj-eonductance in dirty metals 
contain a small parameterlia l/pphmpi an d the absence 
of this parameter in the noise magnitude implies that 
quantum interference is irrelevant to shot noise. Afcthe 
values of higher cumulants obtained by Lee et alfcl ei- 



ther do not contain the quantum-interference prefactor 
as compared to the second cumulant, they can be also 
obtained semiclassically. T_bis conclusion is supported 
by recent numerical results, L3 which were obtained in a 
semiclassical model with exclusion principle. 

The present paper extends the semiclassical Boltz- 
mann - Langevin approach to calculations of the third 
and fourth cumulants of current for diffusive-metal con- 
tacts. We show that the main contribution to these quan- 
tities still results from second-order correlation functions 
of the Langevin sources in the Boltzmann - Langevin 
equation, whereas the effect of higher-order correla- 
tions of these sources is negligibly small in the diffusive 
case. As the result involves several hierarchically coupled 
second-order correlators, this approach may be termed 
" cascade" . 

The Boltzmann-Langevin equation for fluctuations 
reads 
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(1) 



where 51 is the linearized collision integral and SJ ext is 
the random extraneous flux. Fluctuations of physical 
quantities are expressed in terms of the fluctuations of 
the distribution function 5f. For example, the fluctua- 
tion of current density is given by an expression 



Sj(r,t) = e^v<5/(p,r,i). 



(2) 



In principle, the nth-order correlation function of any 
physical quantity may be calculated provided that the 
correlation functions of the extraneous sources 5J ext are 
known up to nth order- 

Kogan and ShulmanB calculated the second-order cor- 
relation function of the extraneous sources based on very 
simple physical considerations. The collision integral in 
the Boltzmann equation is of the form 



J(p, r, t) = J2 [J(P -> P') - J(P' -» P) 



(3) 
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where J(p — > p') and J(p' — ► p) are the outgoing and in- 
coming scattering fluxes from and to state p. It is natural 
to assume that all scattering events are local in space and 
the scattering fluxes between different pairs of states are 
statistically independent. Furthermore, since scattering 
events at different instants of time are independent, the 
scattering between each pair of states presents a Pois- 
sonian process whose second-order correlation function 
(6J(ti)8J(t<i)) is proportional to its average rate 

J(p - p') = W(p, p')/(p, r, t)[l - f(p\ r, t)}, (4) 

where W(p, p') is the classical scattering probability and 
the factor 1 — / takes into account the Pauli principle. 
Hence the correlation function of extraneous sources re- 
latccLwith randomness of scattering may be written in a 
forma 

(SJ ext ( Pl , n, ti)8 J ext ( P 2, r 2 , t 2 )) = 5( ri - v 2 )5{h - t 2 ) 



J PlP2 



]T[ J (Pi "> P') + J (P' -> Pi)] - J (Pi -» Pa) 



p' 



J(P2 -» Pi) 



(5) 



The correlation function is positive if pi = P2 and nega- 
tive if pi 7^ P2- In the latter case the only correlation is 
possible through scattering from p! to p 2 and vice versa, 
but this results in fluctuations of the corresponding oc- 
cupancies of opposite signs. 

In the case of a diffusive metal and low frequencies, the 
Boltzmann - Langevin scheme is greatly simplified. It is 
convenient to introduce a fluctuation. p£ electric potential 
86 and random extraneous currentEj'Ej 



6r t (r,t)=erJ2^J ext (p), 



(6) 



where r is the elastic scattering time. The fluctuation of 
current density is given by an expression 



8} 



-aVScj) + 8j ext , 



(7) 



where a is the metal conductivity. Because pile-up of 
charge is forbidden in the quasi-static limit, one easily 
obtains from the condition of current conservation that 

aV 2 56 = Vj ext . (8) 

The correlation function of extraneous currents 

(tjT^MW^^M)) = 2*8*1,8(1! - r 2 )S(tt - i a ) 



J tfc/(e,n)[l-/(e,n)] 



(9) 



is easily obtained from Eq. (||). Expressions (Q) - (j|) 
form a complete set of equations for calculating the sec- 
ond cumulant of current. 
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FIG. 1. Typical contributions to (a) second, (b) third, 
and (c) fourth cumulants of current in a diffusive metal. 
The arrows correspond to current fluctuations at different in- 
stances of time. Black circles correspond to the correlators of 
Langevin sources, empty triangles and squares, to their first 
and second functional derivatives. 



It is natural to use the same notion of independent 
scattering events for determining the correlation func- 
tions of Langevin sources of higher order. According to 
the ideas of Kogan and Shulman, the Langevin source is 
a fluctuation of the difference between the outgoing and 
incoming scattering fluxes into a given state (p, r). Since 
each scattering flux between a pair of states is assumed 
to be a Poissonian process, all its cumulants are pro- 
portional to the average scattering rate 1/t. Hence the 
third- and fourth-order correlation functions of extrane- 
ous currents (Q) are proportional to r 3 (see Appendix for 
details). As will be shown below, the contribution from 
these correlations is negligibly small as compared to an- 
other mechanism related with second-order correlations 
of Langevin sources. These correlations may be called in- 
direct, or "cascade". The point is that the second-order 
correlation function of extraneous sources is a functional 
of the distribution function / and hence may also fluctu- 
ate as the current itself. One may roughly imagine these 
fluctuations as fluctuations of Nyquist noise of a resistor 
caused by fluctuations of its temperature. A fluctuation 
8f results in a fluctuation 

*<«S rt (ri,*i)«f rt (*2,*2)> = 2a8^S(n -ra)$(ti -t 2 ) 



J tfe[l-2/(e,n)]*/(e,ri,ti). 



Note that the time scale for the fluctuation 8f is set by 
a characteristic relaxation time and hence they are slow 
as compared to the duration of a scattering event. As 8f 
is also correlated with fluctuations of other quantities, 
this "noise of noise" results in higher-order correlations 
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of currents. Hence even Gaussian extraneous sources may 
produce non-Gaussian fluctuations. 

In a symbolic form, the expression for indirect third- 
order correlations may be written in a form 

(Shsi 2 si 3 ) = p 123 | y (w>| , (io) 

where P123 denotes a summation over all inequivalent 
permutations of indices (123) and 5(...)/5f 4 denotes a 
functional derivative with respect to f{e4,,T^,ti). The 
products imply a convolution over the arguments of the 
distribution functions with repeating indices. The irre- 
ducible fourth-order correlation function is presented by 
three groups of terms 



8h5I 2 8h5h :) = Pi 



12:il 



S 2 (ShSI 2 ) 



{Sf 5 SI 3 ){Sf e Sh) 



SjShSh) S(5f 5 5I 3 ) 
Sf 5 6f 6 



S(ShSI 2 ) S(SI 3 Sh) 



Sf 6 



(11) 



The corresponding contributions are schematically illus- 
trated by diagrams in Fig. 1. 

To evaluate expressions ( |To| ) and (fy]), one must know 
the second-order correlation function of fluctuations Sf. 
In the quasi-static limit, the Boltzmann - Langevin equa- 
tion may be easily transformed into a stochastic diffusion 
equation for the fluctuation Sf(e, r, t) 



DV 2 5f{e,v) = V<SF e 



(12) 



where 



5W xt (e) = ^ £ ^J^(p)S(s p - a), (13) 



so one obtains from (n|) that 



D 



(5F™\e, r)SF^(e', r')> = 2^-S(r - r')6(t t>) 
S(e-e')S aP f(e)[l-f(e)}, 



(14) 



(6F^(e, v)8 ] f\v')) = 2eD5{v r')S(t - t') 

aS aP f(e)[l-f{e)], (15) 

where D is the diffusion coefficient and Nf is the Fermi 
density of states. 

Consider the case of a quasi-one-dimensional diffusive 
wire of length L connecting two massive electrodes with 
a voltage drop V across it. In this case, all quantities 
may be considered as depending only on the longitudinal 



coordinate x along the wire. The solution of Eq. (|lj) is 
of the form 



Sf(x) 



DSn 



dx'K{x,x') / d J r±5F* xt (e,x,r ± ), 



K(x, x) = 9{x — x) — x/L, 



(16) 



where So is the area of the cross-section of the wire and 
r± are the transverse coordinates. Because the pile-up 
of charge in the contact is forbidden in the quasi-static 
limit, the fluctuation of current is independent of x and 
is obtained by averaging the extraneous current over the 
contact volumelll] 



si= JtJ d2r±sjr 



(17) 



The low-frequency Fourier transform of the correlation 
function of 61 (ti) and SIfa) is given by 



{5I(u!l)SI(u)2)) — S(ui + L0 2 ) 



x r~l J * dx J de f( £ > x ')[ 1 ~f( e > x> )]> 



(18) 



where R is the resistance of the contact. By giving the 
distribution function a variance 5f(e, x) , one easily ob- 
tains the functional derivative of Eq. (|l|) with respect 
to f(e, x, 0J3) in the form 

5(51(^)61^)) 4tt 

— — — =6(l} 1 +u 2 -u 3 — 1-2/ £, x )}. 

5f(e,x,u)3) RL 

(19) 

The second functional derivative of the same quantity is 
just 



6 2 (5I(lu 1 )5I(u J2 )) 

8f{£3, X 3 ,UJ 3 )5f(e 4 , X 4 , W 4 ) 



-5{<jJ\ + LU 2 — UJ 3 - L) 4 ) 



x 5(e 3 - e 4 )5(x 3 - x 4,)^- 



(20) 



The right-hand side of Eq. ( |10| ) allows three inequivalent 
permutations of indices. Hence the low- frequency Fourier 
transform of the third cumulant of current 

S 3 {cj 1: oj 2 ) = I d(h - t 3 ) ( d{t 2 - t 3 ) 



x exp[iwi(<i - t s ) + ioj 2 (t 2 - * 3 )](W(ii)W(ta)«/(t 3 )> 
is given by an expression 

S 3 (0,0) = l2 J[fj2- J d£ j dx J dx ' \ l - 2 f( £ i x )] K { x , x ') 
xf(e,x')[l-f(e,x')}. (21) 
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The first, second, and the third terms in the right-hand 
side of Eq. ( |TT| ) allow 6, 12, and 3 inequivalent permu- 
tations of indices, respectively. Hence the low-frequency 
Fourier transform of fourth cumulant of current 

S4,(ux,u>2,u 3 ) = / d(ti-t 4 ) d(t 2 -U) d(t 3 -t 4 ) 



x exp[zwi(fi - t 4 ) + zcj 2 (^2 — ti) + iW3(*3 — £4)] 

x(: 6I{ti)SI{t 2 )SI{t z )6I{U) :> 
may be written in a form 

5 4 (0, 0, 0) = 6S4-1 + 12^4-2 + 3S 4 -3, (22) 

where 

2 L L 

S4—1 = — 16 ^^ 3 J de J dx J dx\K{x,xi) 


x/(e,a:i)[l - f{e,x{)] 

L 

x J dx 2 K(x,x 2 )f(e,x 2 )[l- f(s,x 2 )}, (23) 


Si-2 = 8-^3 J de J dx[l - 2f(e,x)] 



x J dxiK(x,x 1 )[l-2f(e,x 1 )} 


L 

x J dx 2 K(x,x 2 )f(e,x 2 )[l- f(s,x 2 )}, (24) 




and 



54-3 = 8-^3 /" de J dx 1 [l-2f{s,x l )] 



5 10 15 20 25 30 35 40 

eV/T 

FIG. 2. Voltage dependence of the fourth cumulant of 
current for a long diffusive contact 

The distribution function inside the contact is given by 
an expressionEJ 

f(e,x)=(l-^)f (e + eV/2) + ^f (e-eV/2), (26) 

where fo(s) is the Fermi distribution function. The in- 
tegration in Eqs. ( pl| ) - ( p5| ) is easily performed. One 
obtains that the third cumulant is equal 



5, = 



1 e eVcosh(eV/T) + 12Tsinh(eV/T) - 13eV 



15 R 



cosh(eVyT) - 1 



(27) 



This is essentially the same expression that was obtained 
by Gutman and GeferJ-3 using the nonlinear tr-model. It 
gives 



1 e 2 V 



15 R 

at high voltages in agreement with Lee et alE and 



5 a = ~ 



le 2 V 



3 R 



(28) 



(29) 



in the low- voltage or high-temperature limit E 

The fourth cumulant of current is given by an expres- 
sion 



S l 



1 e 2 



1 



420 i?sinh J (eV72T) 



eV cosh 



3eV 



x J dx 2 [l - 2/(e, x 2 ) 



x / dxK{x 1 ,x)K(x 2 ,x)f{e,x)[l - f(e,x)]. (25) 
o 



-20Tsinh 



3eV 
2T 



-(M7\inl. [ |^ 



(30) 



It gives 
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1 



in the hig, 
Lee et al 
inEq. (fj2| 



105 R < 31 > 

voltage limit in agreement with the results of 
In equilibrium at V = only the last term 
is nonzero and 



(32) 



2 e 2 T 

This result seems to be new. The positive sign of S4 in 
equilibrium is in a qualitative agreement with the shape 
of distribution function-pf current numerically obtained 
by Belzig and NazarovES The behavior of Si(V) is shown 
in Fig. 2. The fourth cumulant changes its sign at 
eV ~ 20T, hence the distribution of current changes from 
super-Gaussian to sub- Gaussian .as the voltage increases. 

It has been proposed recentlyO that the third cumu- 
lant of noise may be used for determining the effective 
charge of quasiparticles. In other words, it may give the 
same information as the shot noise even at low voltages 
eV <C T. Since the fourth cumulant is nonzero even 
at zero voltage, one may obtain a nontrivial information 
about the system from measurements of equilibrium fluc- 
tuations. Measurements of fluctuations of noise power, 
which are termed "the second spectral density" and are 
closely related with the Fourier transform of the fourth 
cumulant J13 were successfully implemented by a number 
of authors who studied the 1// noise£H~0 

In summary, we have shown that higher cumulants 
of current in diffusive-metal contacts may be described 
by the semiclassical Boltzmann - Langevin theory. The 
semiclassical fluctuations appear to be^on-Gaussian 
both in equilibrium and at high voltages.B3 The present 
approach may be also applied to semiclassical non- 
diffusive systems. In this case, one has to take into 
account also higher-order local correlators of extraneous 
sources. 

The author acknowledges a helpful discussion with Sh. 
M. Kogan. This work was supported by Russian Foun- 
dation for Basic Research, grant 01-02-17220, and by an 
INTAS Open grant, and by an INTAS Open grant 2001- 
1B-14. 

APPENDIX 

The third central moment of 8 J ext is easily found from 
the same considerations as the second moment. Much 



like the second-order correlator, it is local in space. As 
the collision integral involves only single-particle tran- 
sitions between pairs of states, the triple correlator 
((5J ea:t (pi)<5J e:,:t (p2)<5J e2: *(p3)) is nonzero only if at least 
two of these states coincide. We also make use of the 
fact that the third-order correlator of a scattering flux 
from one state to another, which is assumed to be Pois- 
sonian, is proportional to the average scattering rate be- 
tween these states. The sign of the contribution from 
each scattering process should be determined by multi- 
plying the signs of occupancy changes in states pi, P2, 
and P3 when the corresponding scattering event takes 
place. Hence the correlator may be written in a form 

{5 J** ( Pl , n , h ) 5J ext (pa , r 2 , t 2 ) SJ ext (p 3 , r 3 , t 3 )) 



= <J(n - r 2 )S(v 2 - v 3 )6{ti - t 2 )5{t 2 - t 3 ) 



<S P1P2 S P2P3 [J(p' -> Pi) - J(pi p')] 
I p' 



+W2 [J(Pl -> P3) - J(P3 -> Pl) 



+ <W 3 [ J (P2 -> Pl) - J (Pl -> P2)] 



J PlP3 



[J(pi -> p 2 ) - J(p 2 Pi)] 



(Al) 



where the scattering fluxes J are given by Eq. (Q). It 
is easily verified that this correlator conserves the total 
number of electrons at a given point. 

Consider now the fourth-order irreducible correlator of 
extraneous sources. The term "irreducible" means that 
only the part of the correlator that cannot be decoupled 
into products of second-order correlators is considered. 
This correlator is zero for a Gaussian noise but is pro- 
portional to the average rate for a Poissonian process. 
Using the same ideas, the irreducible correlator may be 
written in a form 



5J ext ( Pl , n, t 1 )SJ ext (p 2 , r 2 , t 2 )SJ ext (p 3 , r 3 , t 3 )6J ext (p 4 , r 4 , U) :} 



S(ri - r 2 )6(r 2 - r 3 )<5(r 3 - r^fa - t 2 )5{t 2 - t 3 )S(t 3 - U) 

X S *piP2 (S P2P3 ( 5p3P4 [ J (P' Pl ) + J (Pl P')] - (S P1P2 5 P2P3 [-/(Pl ~> Pi) + J(P4 ~> Pl) 

I P' 



^PiP2^ P 2P4 [J(Pi Pa) + ^(Ps Pi)] - (S PlP3 S P3Pi + <5 P2P3 '5p3P4) [J(Pi P2) + J(P2 Pl)] 
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+ <5piP2<W4 [ J (Pl P3) + ^(P3 -»■ Pi)] 



-(*P! 



P3^P2P4 



^PlP4^P2P3 ) 



[J(pi -> p 2 ) + J(P2 -> Pi)] 



(A2) 



In the diffusive limit, it is convenient to pass from ex- 
traneous fluxes in the Boltzmann - Langevin equation 
to extraneous currents in the drift-diffusion equation by 
means of Eq. @ . The third-order correlator of extrane- 
ous currents is of the form 



<J(ri - r 2 )<5(r 2 - r 3 )S(h - t 2 )6(t 2 - t 3 ) 



pp' 



x[J(p'-p)- J(p-p')]. 



(A3) 



Though the scattering fluxes J are proportional to 1/r, 
the bracketed difference vanishes in equilibrium and is 
proportional to the anisotropic part of the distribution 
function in the case of applied bias. As this part is pro- 
portional to t in the diffusive approximation, Eq. (A3) 
is at least of the order of t 3 . 

The fourth-order irreducible correlator of extraneous 
currents is of the form 



16 
15 



e^Lp^Oi - r 2 )£(r 2 - r 3 )S(r 3 - r 4 ) 



x6(ti - t 2 )S(t 2 - h)8(t 3 - t4)(S a /3S 7 s + SajS^s + 8 a s$i3j) 



de/(e,ri)[l - /(e,ri)], 



(AA) 



where a is the conductivity. Hence the fourth-order cor- 
relator is also of the order r 3 . 



Y. Blanter and M. Biittiker, Phys. Rep. 336, 1 (2000). 
L. S. Levitov and G. B. Lesovik, Pis'ma Zh. Eksp. Teor. 
Fiz. 58, 225 (1993) [JETP Lett. 58, 230 (1993)]. 
' H. Lee, L. S. Levitov. and A. Yu. Yakovets, Phys. Rev. B 
51, 4079 (1995). 

D. B. Gutman and Y. Gefen, |cond-mat/0201007 



5 A. Kamenev, A. Andreev, Phys. Rev. B 60, 2218 (1999). 

6 M. J. M. de Jong, Phys. Rev. B 54, 8144 (1996). 

7 Ya. M. Blanter, H. Schomerus, and C. W. J. Beenakker, 
Physica E 11, 1 (2001). 

8 Ya. M. Blanter and E. V. Sukhomkov, Phys. Rev. Lett. 84, 
1280 (2000). 

9 Sh. M. Kogan and A. Ya. Shul'man, Zh. Eksp. Teor. Fiz. 
56, 862 (1969) [Sov. Phys. JETP 29, 467 (1969)]. 

K. E. Nagaev, Phys. Lett. A 169, 103 (1992). 

1 C. W. J. Beenakker and M. Biittiker, Phys. Rev. B 46, 
1889 (1992). 

2 A. A. Abrikosov, Fundamentals of the Theory of Metals, 
North Holland, Amsterdam, 1988. 

3 P.-E. Roche and B. Doucot, |cond-mat/020214c . 



K. E. Nagaev, Phys. Rev. B 57, 4628 (1998). 
' E. V. Sukhomkov and D. Loss, Phys. Rev. Lett. 80, 4959 

(1998); Phys. Rev. B 59, 13054 (1999). 
' W. Belzig and Yu. V. Nazarov, Phys. Rev. Lett. 87, 067006 

(2001). 

L. S. Levitov and M. Reznikov, :ond-mat/01 11057 . 
' Sh. Kogan, Electronic Noise and Fluctuations in Solids 

(Cambridge University Press, Cambridge, 1996). 
' G. A. Garfunkel, G. B. Alers, and M. B. Weissman, Phys. 

Rev. B 41, 4901 (1990). 

1 N. E. Israeloff, G. B. Alers, and M. B. Weissman, Phys. 
Rev. B 44, 12613 (1991). 

C. E. Parman, N. E. Israeloff, and J. Kakalios, Phys. Rev. 
Lett. 69, 1097 (1992). 

: Both of thesa statements are in contrast with the claim 
made in Ref.u The conclusion of these authors was based 
on the fact that the third cumulant of current cannot be 
obtained correctly with purely local extraneous sources in 
the standard version of the Boltzmann-Langevin approach, 
which does not involve the above cascade mechanism. How- 
ever even this standard version yields non-Gaussian fluctu- 
ations in the quasi-ballistic regime (see Appendix). 



G 



